Suppression of static stripe formation by next-neighbor hopping 
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We show from real-space Hartree-Fock calculations within the extended Hubbard model that next- 
nearest neighbor (t 1 ) hopping processes act to suppress the formation of static charge stripes. This 
result is confirmed by investigating the evolution of charge-inhomogeneous corral and stripe phases 
with increasing t' of both signs. We propose that large t' values in YBCO prevent static stripe 
formation, while anomalously small t' in LSCO provides an additional reason for the appearance of 
static stripes only in these systems. 



The presence of charge-inhomogeneous and striped 
phases as candidate ground states for the cuprates 
continues to be one of the most contentious issues 
in high-temperature superconductivity (for a review 
and references see Ref. [j]). Static stripes have been 
observed experimentally, but to date only in rare- 
earth-doped La2- a; Sr 2 ;Cu04 (LSCO) systems with the 
low-temperature tetragonal (LTT) structural distortion. 
Theoretical explanations for stripes fall into three cate- 
gories: i) they are a true ground state intrinsic to short- 
ranged models of the Cu02 planes, ii) they are a com- 
peting excited state stabilized by anisotropy, or iii) they 
emerge when long-ranged interactions are invoked (gen- 
erally to frustrate phase separation). We have recently 
subscribed to the second viewpoint, motivated by the 
direct experimental connection of lattice structure with 
stripe formation and suppression of superconductivity^ 

However, the interpretation of striped phases as a con- 
sequence of lattice-induced hopping and superexchange 
anisotropy in the Cu02 plane leaves some open ques- 
tions. On the qualitative level, these include the absence 
of static stripes in orthorhombic YBa2Cu307_5 (YBCO) 
systems, while quantitatively they include the value of 
the relative anisotropies required to stabilize a stripe 
phase. A key issue long recognized as a source of ma- 
jor differences between cuprate compounds is the shape 
of the Fermi surface.EI This has been measured by angle- 
resolved photoemisaion (ARPES) for Bi 2 Sr2CaCu2jp 8 +5 
(BSCCO)j3 LSCO,y and Nd 2 _ a; Ce a: Cu04 (NCCO)0 sys- 
tems, and is most easily modeled by extended tight- 
binding band structures .□ The primary influence of an ex- 
tended band structure may be encapsulated in the single 
parameter t' governing the next-nearest-neighbor hop- 
ping on the square lattice. An investigation of t' effects 
is required in context of stripes, and is provided here 
within the real-space Hartree-Fock (RSHF) technique!!! 

The effects of next-neighbor hopping on stripes have 
been considered in .-the extended t-J model by exact di- 
agonalization (ED)El and by density-matrix renormaliza- 
tion group (DMRG) calculations.!] Both sets of authors 
reported a suppression of horizontal (or vertical) stripe 
order with increasing \t'\ of both signs, but detailed fea- 
tures of the results were not fully consistent. The authors 



of Ref. also made a systematic comparison of stripe and 
pairing instabilities, finding a strong anticorrelation. A 
t' opposite in sign to t was included in a more general 
study of the extended (t-t') Hubbard model by RSHF,II 
which confirmed the tendency towards stripe suppres- 
sion. The influence of t' on stripes has been considered 
in the context of Fermi-surface geometry,EI and in terms 
of its effects on stripe filling.Eij Finally, extended hopping, 
integrals were used in an ED study of doped stripes,EJ 
where their role was primarily to avoid phase separation. 
We return below to a more detailed discussion of these 
results. 

The extended, anisotropic Hubbard model is given by 
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We retainjn Eq. (|l|) the possibility 
of anisotropic nearest-neighboru and next-neighbor hop- 
ping, t x / t y and t' x+y ^ t'x-y However, the symmetry 
of next-neighbor terms in this one-band model is main- 
tained in the LTT phase of LSCO, and for most of the 
study to follow. We apply the HF decomposition of the 
Hubbard term in Eq. (lj) and seek self-consistent solu- 
tions for the static charge and spin configurations. 

A detailed analysis of the RSHF technique was pro- 
vided in a_ recent study of anisotropy effects on stripe 
formation!!! The solutions were characterized as functions 
of the "intrinsic" system parameters, by which is meant 
the ratio U/t, the (hole) filling x, the temperature T, and 
the hopping anisotropy t x 7^ t y , and for their dependence 
on the "extrinsic" parameters system size and geometry, 
boundary conditions (BCs), and commensuration of x 
with system size. As a result of this investigation, we 
have chosen two parameter sets representative of inho- 
mogeneous spin and charge configurations whose evolu- 
tion we will follow on varying t' in the extended Hubbard 
model. Hereafter we take the energy unit to be t = 1. 

We focus on 12x12 systems with open BCs, and take 
U/t = 5 to give line-like charge structures. The cases we 
consider are: I - isotropic hopping, t x = t y = 1, with 
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FIG. 1. Ground-state charge distributions for Hubbard 
model on a 12 x 12 cluster with U = 5, t' = 0, and open BCs. 
In (a) x = 1/8 and t x = t y = 1, while in (b) x — 1/6, t x — 0.9 
and t y = 1.1. In these figures the hole density is scaled by 
radius, and the largest circles correspond to (rii) = 0.644, or 
35.6% hole doping of the site. 

hole doping x = 1/8, where the ground state is a corral, 
or closed loop of diagonal, antiphase domain walls; II - 
anisotropic hopping, t x = 0.9, t y = 1.1, with hole doping 
x = 1/6, where the ground state consists of uniform, 
filled, antiphase stripes. The charge configurations with 
t' = are reproduced!! in Figs. 1(a) and (b) respectively 
for cases I and II. The site- or bond-centered nature of 
the charge structures is not important for studying the 
influence of t' . With regard to stability against small 
changes in t', both corral and stripe solutions are found to 
have only minor, quantitative alterations when < 0.05 
for both signs. 

Fig. 2 illustrates the effects of increasing a next- 
neighbor hopping t' with the same sign as t. This results 
in a change of shape of the Fermi surface, which remains 
closed around the T point but expands in the directions 
(k, ±k) while contracting in the (±A;, 0) and (0, ±/c) direc- 
tions. For the corral solution we see at t' = 0.1 [Fig. 2(a)] 
a trend towards breaking of the domain line into small 
clusters, followed before t' = 0.2 [Fig. 2(b)] by a com- 
plete evaporation of the inhomogeneous charge structure. 
The solution here and for higher is a charge-uniform, 
"metallic" phase. For the stripe solution the behavior is 
similar, in that the stable stripe phase is lost in the range 
0.1 < t' < 0.2 in favor of the uniform phase, which then 
persists to high values of t' . The spin configuration (not 
shown) in the uniform phase is found to be commensu- 
rate antiferromagnetism, although with an ordered mo- 
ment suppressed in comparison to the antiferromagnetic 
(AF) regions of the corral and stripe solutions at t' = 0. 

This result has a straightforward interpretation in 
terms of the metallic state maximizing the hopping ki- 
netic energy. Positive t' terms cause direct frustration 
of local charge and spin alignment by competing with 
the t term, resulting in a uniform distribution. One may 
also consider a band picture, where for this sign of t' 
the bandwidth of the lower Hubbard band in the spin- 
density wave (SDW) state at half-filling is increased, sug-» 
gesting that doped hole-like carriers favor delocalization.B 
Our result and interpretation are fully consistent with 
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FIG. 2. Ground-state charge distributions for Hubbard 
model on a 12 x 12 cluster with (7 = 5 and open BCs. In 
(a)-(c) x — 1/8 and t x — t y — 1, while in (d)-(f) x — 1/6, 
t x = 0.9 and t y = 1.1. t' = 0.1 in (a) and (d), t' = 0.2 in (b) 
and (e), and t' = 0.3 in (c) and (f). Charge scale as in Fig. 1. 

DMRG studies of the t-J model,@ which for t' > show 
a systematic loss of stripe charge order. These also find 
an enhancement of hole pairing, which is further con- 
sistent with our result because superconductivity would 
be expected as the leading low-temperature instability 
of the isotropic, metallic state. ED calculations^] in this 
regime differ in that they report enhanced stripe order for 
< t' < 0.2, albeit for diagonal stripes at smaller (phys- 
ical) values of J/t. The authors interpret this result in 
terms of the explicit next-neighbor hopping required to 
cancel the intra-sublattice hopping generated implicitly 
by the t term, and in accord with the present findings 
expect that carrier localization and stripe-formation are 
maximized when the sum of these two contributions to t' 
is minimized. A DMRG studyS which tends to support 
this result remains in minor contradiction to Ref. and 
the resolution is presumably to be found in commensu- 
ration effects on the differing system sizes used. 

Fig. 3 illustrates the effects of increasing \t'\ in the 
physical regime for the cuprates, where t' has the oppo- 
site sign to t. This results in the change from a closed 
to an "open" Fermi surface around the T point. For the 
corral solution we see at t' = —0.1 and —0.2 [Figs. 3(a,b)] 
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FIG. 3. Ground-state charge distributions for Hubbard 
model on a 12 x 12 cluster with U = 5 and open BCs. In 

(a) -(c) x = 1/8 and t x — t y = 1, while in (d)-(f) x = 1/6, 
t x = 0.9 and t y = 1.1. t' = -0.1 in (a) and (d), t' = -0.2 in 

(b) and (e), and t' = —0.3 in (c) and (f). Charge scale as in 
Fig. 1. 

a reinforcement of undisturbed AF order in the center of 
the system, and retention of the diagonal domain walls, 
until a more complex behavior sets in at t' = —0.3 
[Fig. 3(c)]. For the stripe solution there is a definite 
crossover from horizontal stripes at small \t'\ to diagonal 
stripes at t' = —0.1 and —0.2 [Figs. 3(d,e)], again preced- 
ing a less clear structure at t' — —0.3 [Fig. 3(f)]. For both 
parameter sets, the high-i' configurations shown in Fig. 4 
are rather complex in shape, only weakly inhomogeneous 
in charge structure and weakly but commensurately AF 
in spin structure. 

For t' Jt < 0, the narrowing of the lower SDW Hubbard 
band suggests reduced carrier mobility and thus an en- 
hanced tendency toward charge inhomogeneity. In a local 
hopping picture, all next-neighbor processes cost kinetic 
energy, and so are suppressed!! by i) aligning spins fer- 
romagnetically on diagonal bonds, promoting AF order, 
and ii) aligning holes on diagonal bonds, promoting diag- 
onal domain walls. Both of these tendencies are clear in 
Fig. 3. At higher t' this trend cannot be sustained, as the 
kinetic energy cost of failure to delocalize is too great, and 
for —0.4 < t' < —0.3 the kinetic energy associated with 




(a) (b) 
FIG. 4. Ground-state charge distributions for Hubbard 
model on a 12 x 12 cluster with U = 5, t' = 0.6, and open 
BCs. In (a) x — 1/8 and t x — t y — 1, while in (b) x = 1/6, 
t x — 0.9 and t y = 1.1. Charge scale as in Fig. 1. 

t' becomes negative despite the sign, as shown in Fig. 5. 
This is possible by a change in the coefficients of the HF 
wave function such that on average the sign of the diag- 
onal overlap is negative. While such an alteration is not 
readily visualized in terms of charge configurations, its 
consequences are the weakly inhomogeneous structures 
in Fig. 4. The tendency to carrier delocalization also for 
this sign of t' is in accord with the results of Ref. [To] for 
the hole filling of an isolated ^(gipe. 

In the DMRG calculations J3i3 increasing \t'\ up to 0.3 
introduces a doubling of the charge periodicity From the 
above considerations we may offer a change from horizon- 
tal to diagonal stripes as a possible consistent explanation 
of this result, given the small numbgr of holes involved. 
The observed suppression of pairingu is certainly consis- 
tent with the expectation of charge localization accom- 
panying the formation of static^ inhomogeneous charge 
structures. In the ED results, □ both vertical and di- 
agonal stripe correlations are strongly suppressed by t', 
whereas pairing correlations remain appreciable. Other 
RSHF studies of the extended Hubbard modeH confirm 
a significantly reduced tendency to stripe formation in 
this parameter regime. 

All of the above results are for hole-doped systems. In 
the electron-doped regime, we find by RSHF an exact re- 
production of the charge configurations in Figs. 1-4, with 
the important proviso that the sign oit' be reversed. This 
symmetry emerges from the band and local physics argu- 
ments invoked above, and can be seen in Fermi-surface 
shapes, but most importantly can be proven rigorously 
from the invariance of the Hamiltonian [Eq. (1)] under 
particle-hole transformation combined with sign reversal 
c,; — > — a on one sublattice. 

In summary, next-neighbor hopping terms t' lead to 
a rapid suppression of the t' = inhomogeneous charge 
structures in the Hubbard model. For hole- (electron-) 
doped systems the primary influence of positive (nega- 
tive) t' (Fig. 2) on stripes lies in its delocalization effect, 
while that of negative (positive) t' (Figs. 3,4) is destruc- 
tion of the spin registry and antiphase domain wall na- 
ture. This latter result suggests the impor.tance of AF 
order in stabilizing static horizontal stripes .E_2l 

In cuprates (t'/t < 0), the effective t' is intrinsically 
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FIG. 5. Average kinetic energy -Ekin per diagonal bond as 
a function of next-neighbor hopping t' opposite in sign to t. 

smaller for LSCO and related monolayer systems than for 
YBCO, BSCCO and the related mono-, bi-, and trilayer 
systems. This result may be obtained by simple fitting of 
tight-banding models to the Fermi surfaces measured by 
ARPESB and is also confirmed by reduction (downfold- 
ing) of band structure calculations.E3 In addition to the 
number of layers, the apical oxygen atoms are seen to be 
an important factor contributing to the difference. Be- 
cause deviations of the Cu-O-Cu bond angle from 180° 
reduce orbital overlap, the LTT and low-temperature or- 
thorhombic (LTO) distortions of the LSCO structure also 
act to reduce \t'\ in comparison with unbuckled cuprate 
layers. At the tight-binding level, when t and t' are the 
only parameters in the extended hopping band structure, 
t' values of order —0.5t are required to model the open, 
YBCO-like Fermi surface, while the closed LSCO surface 
is reproduced by values of \t'\ < 0.2. 

We have considered previously the possible key role 
of the LTT distortion, which is the ground-state lattice 
structure only in Nd- and Eu-doped LSCO compounds.EJ 
Only the LTT distortion produces anisotropic t terms, 
and to the best of our knowledge there remain no reports 
of static stripes in any non-LTT LSCO structures, or 
in any other cuprates. A natural question then arises 
concerning YBCO, where the presence of 6-axis chains 
leads to a planar structural anisotropy. Our results offer 
an additional explanation for the absence of static stripes 
in YBCO: the large effective t' values [Fig. 4] suppress the 
formation of charge-inhomogeneous structures which are 
candidate ground states at smaller t' . Further, t' may be 
expected to be an important factor contributing to the 
absence of static stripes in the majority of hole-doped 
cuprate materials, where the Fermi surfaces of the Cu02 
planes are open. For the electron-doped systems, we have 
shown these to be in the regime where stripe formation 
is suppressed by any value of t' , which for these dopings 
favors a uniform, metallic phase. 

We comment briefly that the LTO structure of the 
Cu02 planes in LSCO and BSCCO has an asymmetrical 
t'. In RSHF, and for t' opposite in sign to t, a sufficiently 
large anisotropy t' x+y > t' x _ y yields diagonal stripes in 
the direction x + y. However, the required anisotropics 



are physically unjustified, and no effect of this nature is 
expected in the real material. 

In conclusion, we have studied the effects of a next- 
neighbor hopping term within the extended Hubbard 
model on the formation of static stripes in cuprate sys- 
tems. In a RSHF treatment, t' terms of both signs cause 
a rapid suppression of horizontal stripe formation. When 
t' has the physical sign for cuprates (opposite to t), we 
find for hole-doped systems a window in which diagonal 
stripes are favored, and then a broad regime at large 
t' where the compromise charge configuration is only 
weakly inhomogeneous. For the opposite sign of t', a ho- 
mogeneous, metallic phase is preferred. These properties 
are exactly reversed in electron-doped systems. 

In contradiction to observations of superconductivity, 
the "physical" case is dominated by charge localization 
and suppression of hole pairing, whereas competition of 
inhomogeneous charge structures with the superconduct- 
ing instability is reduced for t'/t > 0. However, these 
results are consistent with previous studies of extended 
Hubbard and t-J models. We propose that the weak t' 
terms characteristic of the LSCO system in comparison 
with other high-T c superconductors is an important con- 
tributing factor to the appearance of static stripes only 
in these materials, and that stripe formation in electron- 
doped cuprates is suppressed for all t'. 
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